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ABSTRACT 


This paper presents a review of recent work at NASA-LaRC to compare the predictions of transverse 
fuel injector flow fields and mixing performance with experimental results. Various "cold" (non- 
reactive) mixing studies have been selected for code calibration; these include the effects of boundary 
layer thickness and injection angle for sonic hydrogen injection into supersonic air. Angled injection of 
helium is also included. This study was performed using both the three-dimensional elliptic and the 
parabolized Navier-Stokes (PNS) versions of SPARK. Axial solution planes were passed from PNS to 
elliptic and elliptic to PNS in order to efficiently generate solutions. The PNS version is used both 
upstream and far downstream of the injector where the flow can be considered parabolic in nature. The 
comparisons are used to identify experimental deficiencies and computational procedures to improve 
agreement. 
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area enclosed by the a = 0 curve 

injector diameter 

mass flow rate of injected H 2 

mass flow rate of mixed H 2 

total H 2 mass flow rate from flow field integration 

Mach number 

pressure 

total pressure 

penetration to outer edge of mixing region (height of a Ii2 = .00035 or Xj, e = .005 
contour) 

jet-to-air dynamic pressure ratio 

Reynolds number 

turbulent Schmidt number 
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total temperature 

plate temperature 
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y or Y lateral coordinate 

z or Z coordinate normal to plate 

Z m height of oc max contour 

a mass fraction 

otrnax maximum mass fraction 

a R mass fraction defined by Eqn. 1 

a s stochiometric hydrogen mass fraction (0.0285) 

P air mass fraction (1-a) 

8 boundary layer thickness 

rj OT mixing efficiency 

x\^ m mixing efficiency with a s = 0.01425 

(r\ln)adj mixing efficiency with a s = 0.01425 and assuming missing fuel to be unmixed 

p density 

9 in j injection angle 


INTRODUCTION 

Efficient hypersonic flight at speeds above Mach 6 require the use of supersonic combus- 
tion ramjet or scramjet engines. Research efforts to establish a technology base for scramjet 
engine components have been underway at NASA Langley Research Center (LaRC) since the 
mid 1960’s. This research was focused on the Mach 4 to 7 speed range, the range within the 
enthalpy limits of combustion heated or arc heated facilities, and directed toward the inlet and 
combustor components of the engine. Combustor technology development [Ref. 1] was 
directed at understanding the requirements for and performance of fuel injection, mixing and 
reaction (ignition and flame holding) in the scramjet combustor environment, and incorporated 
a dual approach of experimental and computational technology development. 

Experimental combustor technology development started with single[2-4] and multiple[5- 
7] fuel injector studies in cold flow to characterize fuel mixing, followed by reacting coaxial 
flow configurations[8,9] to characterize the supersonic combustion phenomena. Direct connect 
combustor experiments[10-16] were used to verify empirical models[17] for the mixing limited 
supersonic combustion process developed from the early mixing studies. Finally, the empiri- 
cal combustor models were integrated into a complete subscale engine design and tested[ 1 8] 
over the Mach 4 to 7 speed range in two engine test facilities. 

Computational fluid dynamics (CFD) technology development started with codes which 
contained the physical modeling necessary to predict the supersonic, chemically reacting flow 
encountered in a scramjet combustor[19]. Because of computer memory limitations, these 
codes were limited to two-dimensional or axisymmetric cases until the mid 1980’s. The two- 
dimensional codes were calibrated to appropriate data for non-reacting fuel injection and mix- 
ing[20], and reacting fuel injection and mixing[21,22] in ducted flows. Various finite-rate 
chemistry models were developed and applied for hydrogen fuel[ 19,23] and silane-hydrogen- 
air mixtures[24]. Recent code applications have evaluated non-reactive three-dimensional fuel 
mixing[25] and reactive fuel mixing in confined simulated combustor flows[25-29]. Even 
though they have provided significant insight into these highly complex flow fields, these stu- 
dies have not provided the essential features of code calibration for scramjet fuel mixing. 
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The objective of this paper is to present comparisons of numerical predictions with 
experimental measurements for non-reacting fuel injection into a supersonic freestream and to 
document efforts to improve the agreement. This study will look at a large data set with the 
intent of identifying experimental limitations and correct procedures to use when exercising 
the codes. The primary index of agreement is overall mixing efficiency of the fuel or simu- 
lated fuel. In the absence of experimental mixing efficiency, the index defaults to fuel mass 
or mole fraction distributions. Secondary values of interest include fuel penetration and 
spreading, wall pressure, and any other measured data available. 

A brief description of the experimental tests chosen for this paper is followed by a brief 
description of the SPARK codes. The computational method used to solve the injection prob- 
lems is discussed along with convergence criteria for the elliptic region. The grid and boun- 
dary conditions are also given. Results from the calculations are then compared with data in 
order to demonstrate the accuracy of the numerical modeling. 

EXPERIMENTAL CONFIGURATION 

Details of the flat plate non-reacting injection cases modeled are provided in Table 1. 
Variables in the table include the jet-to-air dynamic pressure ratio, q ; the freestream Mach 
number, M air \ freestream Reynolds number, Re; downstream injection angle, 0 (n ■ , measured 
from the plate surface; the boundary-layer thickness upstream of the injector normalized by jet 
diameter, 8/D; and the jet diameter in centimeters, D. 


Table 1. Summary of Conditions 


Case 

Injectant 

Q 

M air 

Re(m x ) 

®inj 

5/D 

D(cm) 

Ref. 

1 

hydrogen 

1.0 

4.05 

8.3x10* 

90° 

1.25 

0.2540 

4 

2 

hydrogen 

1.0 

4.05 

6.19xl0 7 

90° 

2.51 

0.1016 

3 

3 

hydrogen 

1.0 

4.05 

8.0x10 s 

90° 

6.25 

0.0508 

4 

4 

hydrogen 

1.0 

4.05 

6.19xl0 7 

30° 

3.40 

0.1016 

7 

5 

helium 

0.6 

3.00 

5.2xl0 7 

30° 

2.20 

0.3175 

30 

6 

helium 

3.1 

3.00 

5.2xl0 7 

30° 

2.20 

0.3175 

30 

7 

helium 

0.27 

3.00 

5.2xl0 7 

15° 

2.20 

0.3175 

30 


Table 2. Experimental Flow Properties 


Air Injectant 


Case 

Pt (MPa) 

1 1 (K) 

Mach 

Gas 

Pt (MPa) 

T t { K) 

Mach 

1 

1.72 

295 

4.05 

H 2 

0.33 

289 

1.0 

2 

1.38 

294 

4.05 

H 2 

0.27 

294 

1.0 

3 

1.72 

303 

4.05 

H 2 

0.32 

298 

1.0 

4 

1.38 

300 

4.05 

H 2 

0.27 

295 

1.0 

5 

.655 

290 

3.00 

He 

0.17 

295 

1.0 

6 

.655 

290 

3.00 

He 

0.86 

295 

1.0 

7 

.655 

290 

3.00 

He 

0.076 

295 

1.0 


Cases 1-4 are from tests performed on a flat plate model (Fig. 1) in a Mach 4.05, 23x23 
cm blow down tunnel. Table 2 presents nominal test conditions for the cold air and sonic 
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hydrogen injectors. The facility was operated at stagnation pressures of either 1.38 or 1.72 
MPa, but the four cases selected all operated with unity jet-to-freestream dynamic pressure 
ratio (see Table 1). Instrumentation was provided to measure the stagnation pressures and 
temperatures of hydrogen and air, plate-surface static pressures and hydrogen concentrations. 
In addition, the H 2 mass flow rate was measured using a sharp-edge orifice meter. Plate 
static -pressure orifices of 0.051 cm diameter were located on the longitudinal center line from 
30D ahead of the injector station to 150D downstream for Case 4. In-stream measurements 
were made using a 0.038 cm ID pitot and a 0.102 cm OD conventional static probe. Gas 
samples were taken from the mixing region by the pitot probe and were analyzed with an on- 
line process gas chromatograph having a cycle time of 60 seconds. 

Historically, mixing efficiency, T| m , is defined as that fraction of the least available reac- 
tant (i.e. 0 2 or fuel) which would react if the fuel-air mixture were brought to chemical 
equilibrium without additional local or global mixing. Thus in fuel rich regions, all of the 
local oxygen is considered "mixed," while in fuel lean regions all of the fuel is mixed. Two 
definitions of mixing efficiency are required— one for flows which are globally fuel rich, and 
one for flows which are lean. For fuel lean flows, as those reported here, 

a R p udA 

r (1) 

m H 2 , total I apudA 

a=0 


m = 


2, mix 




where: a is hydrogen (fuel) mass fraction 


a 


where a < a s 


a R =J 


( 


1-a 

1 -a. 


)oc s 


where a > a s 


. v. 

a s . is H 2 stochiometric mass fraction (0.0285) 

A (y ^ 0 is the area enclosed by zero H 2 contour defining the extent of the mixing region 
rh-H mix is tnixed H 2 mass flow rate 

1 2 , total is tota l H 2 mass flow rate from flow field integration. 

When a max < a s , r\ m equals 1.0. 


Experimental mixing efficiencies for the data from Ref. 3, 4 and 7 are determined by 
integration of the hydrogen and airflow contours as illustrated in Fig. 2. Overlaying the sto- 
chiometric fuel-air contour {a s = 0.0285) on the H 2 and air flow contours to provide integra- 
tion limits, the total hydrogen mixed and mixing efficiency are determined by the integration 
procedure of Eqn. 1. Total hydrogen mixed in the fuel lean region is determined by integration 
of hydrogen flow rate in that region (Fig. 2b), while hydrogen mixed in the fuel rich region is 
determined by integration of the air flow rate within the fuel rich region (Fig. 2c). 


As discussed in Ref. 3, one indication of the overall accuracy of the profile data is com- 
parison of the integrated hydrogen mass flow rate with the metered hydrogen mass flow sup- 
plied to the jet. Figure 3 illustrates this comparison for all data in Ref. 3. Agreement of the 
integrated and metered mass flows improves as X/D increases and q decreases. This improve- 
ment is probably due to the smaller gradients in the concentration and velocity which are asso- 
ciated with the local turbulence level at the downstream stations and lower values of q. It has 
also been suggested that selective sampling by the subcritical pitot probe produces the low 
indicated hydrogen composition[30]. 
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Figure 3. Results of Fuel Mass Flow Contour Integration [Ref. 3] 




Cases 5-7 represent low-angled helium injection experiments performed in a Mach 3, 
11x23 cm blow down tunnel at Virginia Teeh[30] with injection from the lower tunnel wall 
Boundary layer thickness at the fuel injector is about two jet diameters. Geometric details of 
the 15° injector presented in Fig. 4 illustrate the converging and short 1-d constant diameter 
section used to provide the sonic helium jet. Nominal test conditions for the cold air and sonic 
helium injector are presented in Table 2. The facility is operated at a total pressure and tem- 
perature of 0.655+.02 MPa and 290 K respectively and helium is injected either at matched or 
underexpanded pressure. Instream measurements of fuel concentration, flow temperature, pitot 
and static pressure were obtained at axial stations of 20, 40 and 90 jet diameters downstream 
of the injector on the jet centerline. Additional surveys of fuel concentration were obtained at 
off jet centerline locations. The gas sampling probe[31] utilizes a very small internal expan- 
sion tip with an internal diameter of 0.028 cm (one tenth the injector diameter), which has 
sufficient internal expansion and flow suction to swallow the tip shock wave. 

DESCRIPTION OF CODES 

The LaRC SPARK family of CFD codes was used in this study. The two-dimensional 
Navier-Stokes code was originally developed and extensively validated for chemically reacting 
and mixing flows by Drummond[32] and has been expanded to three dimensions and validated 
by Carpenter[33]; this three-dimensional version has additionally been developed into a para- 
bolized Navier-Stokes (PNS) version by Kamath[34], The elliptic SPARK code can use either 
a temporally explicit second-order accurate MacCormack-based finite-difference technique to 
solve the mass, momentum, and energy conservation equations or a fourth-order compact spa- 
tial scheme[33] which provides high spatial accuracy and a convenient method for marching to 
the steady state. In addition, the code has the fourth-order time accurate Gottlieb scheme[33]. 
For mixing results, the source terms in the species continuity equations are set to zero. 

At the time of this work, the SPARK codes contained only the Baldwin-Lomax algebraic 
turbulence model. The formulation used contained the constants as presented in the original 
paper[35] for a thin layer on a flat plate. Freestream turbulence in the tunnel was ignored. 
For Cases 5-7, the turbulent viscosity was limited to 1000 times the laminar viscosity. It was 
found in this investigation that the use of Sc t < 1.0 enabled better agreement between numeri- 
cal results and experimental data. There is both experimental and computational justification 
for the use of turbulent Schmidt numbers less than unity for flows such as that examined in 
this work[36,37]. 


SOLUTION METHODOLOGY 

Due to limited computational resources, the solutions were perfonned in a piece-wise 
elliptic/PNS fashion. The flowfield upstream of the elliptic injection region is modeled using 
the SPARKPNS code in 2-d mode. The injection nearfield region, which is highly three- 
dimensional and characterized by such features as large-scale vorticity, separation ahead of 
and behind the injector, strong shocks and expansions, was computed with the 3-d elliptic ver- 
sion of the SPARK code. The injection region extended from six to ten jet diameters ahead 
of the injector to ten to forty jet diameters downstream of the injector, depending on the case. 
The 3-d SPARKPNS code was used in the farfield with the outflow plane from the elliptic 
code passed as a fixed inflow into the PNS code. A small inconsistency is inherent in this 
strategy which is illustrated in Fig. 5, since the character of the governing equations is 
changed discontinuously from the full Navier-Stokes (FNS) to that of the parabolized approxi- 
mation. However, this was considered an acceptable compromise to achieve computational 
efficiency, and the inconsistency involved in the passing of a plane from the full Navier- 
Stokes elliptic domain to the parabolized version of these equations are believed to be small in 
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Figure 5. Solution Methodology 
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relation to turbulence and grid issues. 

The computations presented here were performed on the Cray supercomputers of the 
Numerical Aerodynamic Simulation (NAS) at NASA Ames and the Cray 2 at LaRC. The 
mixing studies each took about 2-10 hours of CPU time. Convergence as defined in the next 
section occurred at about 10000 iterations. The parabolized computations in this study 
comprised about 15% of the CPU time used and a fraction of the memory. 

CONVERGENCE REQUIREMENTS 

Due to the swirling nature of the flow around and behind the injectors, the flow is at best 
quasi-steady, and there is some small fluctuation of pressure in the elliptic region as well as 
for flow parameters in the outflow plane passed to the parabolized code. Claimed conver- 
gence for an elliptic solution of this nature must be limited in its scope to the unchanging 
mean or averaged values of the parameters of interest such as density ((Ap) max < 5%) and the 
conservation (within 10%) of total fuel for all axial planes. 

GRID AND BOUNDARY CONDITIONS 

Grids for each solution discussed herein are summarized in Table 3. All cases used a 
rectangular grid for both the parabolic and elliptic regions with the downstream parabolic 
region having the same cross-sectional grid as in the elliptic. All meshes were clustered near 
the injector and near the plate. Grid domain represents the number of nodes in the axial, 
lateral and vertical directions. The last 5 columns in the table deal with grid dimensions in 
terms of jet diameter with the origin at the injector center (Fig. 5). 

All inflow boundary values were fixed. No-slip boundary conditions were applied to the 
flat plate, which was either at a constant temperature or adiabatic. All of the circular injectors 
were modeled as rectangles on the mesh lower boundary (the plate). The nodes within the 
orifice had their properties set to the jet properties calculated from the stagnation properties 
listed in Table 2. The computational injector was modeled to match the injector area. The 
computational area and/or the velocity on the jet edge nodes were then reduced to match the 
metered injectant mass flow, thereby simulating a discharge coefficient due to losses near the 
injector walls. All other boundary planes were set as outflow planes, except for those of Case 
4 which had multiple injectors and thus two symmetry planes. (Table 3 also lists the number 
of nodes used to model an injector. At each such node, the jet conditions were applied.) 

Table 3. Grid and Solution Details 


Case 

Domain 

Injector 

Sc t 

T 

A vv 

A/D 

B/D 

C/D 

E/D 

W/D 

1 

41x61x41 

25 

Mill i ir mm 

IE i 

-6 

10 

120 

16 

32 


41x51x51 

20 

0.2 

1 

-6 

10 

120 

16 

16 

2 

51x31x41 

25 



-16 

33 

120 

49 

49 


51x31x21 

25 


Bin 

-16 

33 

120 

49 

49 


41x25x21 

9 


' 

-16 

33 

120 

49 

49 

3 

41x61x41 

25 

0.2 


-6 

10 

120 

16 

32 

4 

41x31x41 

25 

0.2 

1 . 

-6 

10 

120 

14 

6.25 

5 

51x41x41 

35 


adia. 

-6 

10 

90 

16 

16 

6 

51x41x41 

35 


adia. 

-6 

10 

90 

16 

16 

7 

51x41x41 

55 

0.5 

adia. 

-6 

10 

90 

16 

16 
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RESULTS 


This section presents computational results and comparisons with experimental data for the 
seven cases. The computational analysis was performed with several values of turbulent 
Schmidt number, and one case was evaluated using several mesh densities. Details of these 
variations are highlighted. The following section will discuss the implications of these findings 
and give recommendations on applying the SPARK or equivalent CFD codes to fuel injector 
problems. 

Case 1 

Computational results for Case 1 are presented in Figs. 6-9. Two different solutions 
were obtained: one with Sc, = 0.2 and one with Sc, = 0.5. Hydrogen concentration contours 
in the vertical X-Z plane through the jet centerline are given in Fig. 6 for the solution with 
Sc, =0.2. These results illustrate fuel feeding upstream into a region of boundary layer separa- 
tion, the penetration and dispersion of the hydrogen core in the nearfield, and penetration of 
the outer edge of the fuel. Note that the a = 0.10 contour (3.5 times stochiometric) extends 
about 3 jet diameters upstream of the injector centerline, and the .005 contour extends to about 
4 jet diameters upstream. The hydrogen core (i.e. the region of highest concentration extend- 
ing downstream from the initial jet, indicated by the dashed line in Fig. 6a) penetrates to about 
1.2 jet diameters above the flat plate and the 0.40 contour (an arbitrarily selected indicator of 
core breakup) persists about 6.5D downstream. The peak hydrogen concentration moves back 
toward the plate and lies along the surface for axial stations past 10 diameters. Fig. 6b shows 
combined FNS and PNS results and indicates complete mixing by 80D downstream (a<a v 
everywhere). 

Computational and experimental results are compared in Fig. 7 by lateral Y-Z plane 
hydrogen mass fraction contours at the 120D station. The outer edge of the fueled region is 
represented here by the 0.003 H 2 mass fraction contour. The computational results have 
greater penetration P/D (see Fig. 2a) and lateral spreading than experimental results with the 
spreading near the plate W b fD nearly double that of the experiment. However, the computed 
0.010 hydrogen mass fraction contour illustrates slightly (=20%) lower penetration and slightly 
(=25%) greater spreading than the experimental contour. 

Computed vertical jet centerline profiles (Sc, = 0.2) at the 120 jet diameter station are 
presented in Fig. 8 and compared with measured H 2 mass fraction, static pressure, static tem- 
perature and velocity. The profiles are in reasonable agreement, except for the location of the 
peak hydrogen concentration. The dashed curve represents results from a solution with Sc, - 
0.5 in both elliptic and parabolized solutions. Again, the peak concentration is adjacent to the 
wall, but the peak value is nearly double that measured because of reduced mass transport. 
However, penetration to the outer edge of the mixing region is not significantly affected by 
the value of Sc, 

Longitudinal variation in predicted fuel mixing efficiency is presented in Fig. 9 for both 
solutions discussed above. Mixing efficiency presented in this figure is NOT as defined in 
Eqn. 1. because of an error in the CFD code used for this and the other hydrogen cases. Mix- 
ing efficiency, r|4, for Fig. 9 is defined by Eqn. 1, but with 

a s = 0.01425 (2) 
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a) H 2 contours, elliptic region. 


o Experimental 



Figure 6. Longitudinal Hydrogen Mass Fraction Contours; Case 1, Sc t = 0.2. 
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Figure 8. Vertical Centerline Profiles at X/D =120; Case 1 , Sc t = 0.2. 



Figure 9. Fuel Mixing Efficiency Distribution; Case 1. 
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Experimental X[' m determined by integration of Ref. 4 contours at X/D = 120 is presented by 
unflagged symbol. Although this is not the standard definition, the comparison of experimen- 
tal and computational results using this method provides a reasonable percent error in model- 
ing mixing efficiency. For this case, between 10% to 20% of the injected fuel was not 
accounted for in the experimental downstream survey. The flagged symbol represents the 
mixing efficiency if all of the missing fuel is assumed to be unmixed; i.e.: 


( y\m)adj ~ ^ 


a R pud A 
m H 2 ,inj 


( 3 ) 


where a s is defined by Eqn. 2. Computational results are in good agreement with experimen- 
tal results when performed using Sc, = 0.2 and with the adjusted experimental result when 
using Sc, = 0.5. At the 120 diameter station, the calculated mixing efficiency for Sc, = 0.5 is 
only 20% lower than for Sc, = 0.2, whereas a max is more than double for the Sc, = 0.2 case. 
The use of an adjusted mixing efficiency results from concern about quality of this older data 
in light of good agreement with the newer He injection data for Cases 5-7. This adjusted 
mixing efficiency represents the lowest bound of experimental scatter conceivable to account 
for limitations. 

Case 2 


Ten computational solutions were obtained for Case 2 for three different values of the 
turbulent Schmidt number on three different density grids plus one laminar fine grid solution. 
Figures 10-13 illustrate flow details obtained with a fine grid and Sc, of either 0.2 or 0.5. 
Figures 14-15 illustrate the sensitivity of mixing efficiency to various grids and Sc, plus one 
laminar solution. 

Hydrogen contours (Sc, = 0.5) presented in Figs. 10a, b illustrate fuel penetration, decay 
of the core and details of the upstream flow separation. The 0.10 contour (3 times sto- 
chiometric) extends less than 1.5 jet diameters upstream of the jet centerline, and the .005 
contour extends about 3.5 diameters upstream. The H 2 core penetrates about 2D, and the 0.40 
contour persists about 6.2D downstream. Peak concentration moves back to the wall at about 
X/D = 25 (not shown). Comparison with experimental penetration of the 0.010 hydrogen 
mass fraction contour are within 10% at X/D = 7, 30, 60 and 120 (Fig. 10b). 

Hydrogen contours (Sc, = 0.5) presented in Fig. 10c illustrate fuel penetration, decay of 
the core and details of the upstream flow separation for the same jet with the same inflow 
plane including boundary layer profile, but with no turbulence in the elliptic region. As 
expected for a disturbed laminar flow, the upstream separation region has become quite exten- 
sive (It should be noted that this solution is not completely converged). The H 2 core 
penetrates about 2.3D, and the 0.40 contour persists about 8D downstream. Both of these 
values are slightly greater than for the turbulent solution. 

Computed lateral spreading of the fuel at X/D at 120D is illustrated in Fig. 1 1 for each 
value of Sc, for the fine mesh. At Z/D = 4.0, lateral spreading W/D (see Fig. 4) of all com- 
putational contours is greater than experimental values, and higher Sc, produces less spread- 
ing. At Z/D = 0, lateral spreading W b /D of all computational contours is in reasonable agree- 
ment with the experimental results. Fuel penetration (P/D) increases for the lowest value of 
Sc , , but all cases exceed experimental results. 
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Figure 11. Lateral Hydrogen Distribution at X/D = 120; Case 2. 
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Figure 12. Vertical Centerline Profiles at X/D =120; Case 2, Sc, - 0.2. 



Figure 13. Longitudinal Centerline Wall Pressure; Case 2, Sc, = 0.5. 
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Figure 14. Fuel Mixing Efficiency Distribution; Case 2. 
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Figure 15. Fuel Mixing Efficiency at X/D = 120; Case 2. 
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Computed vertical jet centerline profiles ( Sc, = 0.2) at the 120D station are presented in 
Fig. 12 and compared with measured hydrogen mass fraction, static pressure, static tempera- 
ture and velocity. Computed profiles are in excellent agreement with experimental results 
except for the fuel concentration. The magnitude of the peak computed hydrogen concentra- 
tion is in excellent agreement with data but is offset to the wall. The dashed curve illustrates 
the computed H 2 mass fraction using Sc, = 0.5. Peak concentration is about 80% higher than 
for the Sc, ~ 0.2 solution. 

Longitudinal jet centerline wall pressure distribution (Sc, - 0.5) in the nearfield is illus- 
trated in Fig. 13. These results show the extent of jet disturbance on upstream flow, recovery 
pressure in the upstream recirculation region, and low pressure in the wake downstream of the 
injector. Note that the pressure rise extends to 5.5D ahead of the injector, which is 
significantly greater than the "zero" hydrogen contour (a = .005) presented in Fig. 10a. 
Downstream pressure is in reasonable agreement with wall pressure data scaled from experi- 
ments in Ref. 2. Calculated results presented are closer to the injector than could be measured 
experimentally. 

Calculated longitudinal distributions of mixing efficiency for Case 2 are illustrated in Fig. 
14 for fine mesh solutions using Sc, = 0.2, 0.5 and 1.0; and in Fig. 15 for all solutions at X/D 
= 120. Mixing efficiency presented in Fig. 14 is as defined by Eqn. 1; whereas the mixing 
efficiency presented in Fig. 15 was calculated using the incorrect value of stochiometric mass 
fraction presented by Eqn. 2. Experimental data presented on each figure have been adjusted 
to the appropriate definition. Fig. 14 illustrates the effect of turbulence on the mixing 
efficiency distribution. A laminar solution (Sc, = 0.5) in the elliptic region illustrates that tur- 
bulence modeling is not an issue for the first 4 jet diameters downstream of the jet, as both 
the laminar and turbulent solutions have identical mixing efficiencies. However, Sc, has a 
large impact in the downstream mixing, because it is used to determine the turbulent mass 
diffusivity from the turbulent viscosity. For Sc, = 0.2, fuel mixing is complete (r) m = 1.0) at 
about X/D = 80, whereas experimental measurements indicate r\ m reaches unity at about X/D 
= 100. Using Eqn. 3 to adjust mixing to account for all hydrogen injected produces flagged 
data points which are in better agreement with the Sc, = 0.5 computational results. Between 
X/D = 5 and 30, turbulent mixing is about 50% higher than laminar mixing. 

Figure 15 illustrates the effect of mesh density and turbulent Schmidt number on the fuel 
mixing efficiency at X/D = 120. Experimental points were obtained using Eqn. 1 and 2 to 
integrate data from Ref. 3. Adjusted data is from Eqn. 2 and 3 using data from Fig. 3 to 
determine the ratio of integrated total fuel to metered fuel flow rate. Fig. 15 illustrates the 
large impact of reduced turbulent Schmidt number on r\ l m for fine grid solutions. Mixing 
efficiency is seen to increase from about 65% with Sc, = 1.0 to 92% with Sc, = 0.2. Coarse 
grids do not have as large a sensitivity to turbulent Schmidt number. Changing vertical mesh 
from 41 to 21 results in increased predicted mixing for Sc, =1.0 from 65% to 75% and no 
change for Sc, = 0.2. The 41x25x21 mesh is coaser both in the x and y directions, as well as 
having the injector defined by 9 nodes versus 25 for the other cases, and the solutions are 
even less sensitive to Sc, than other grids. Note that all solutions with Sc, = 0.5 predict v\ m = 
0.79±0.04. Note also that the data is in the best agreement with fine grid solutions with Sc, = 
0.2 for this case. 

Case 3 

Only one computational solution was performed for the smallest normal jet, and results 
from that case are illustrated in Figures 16-18. The solution utilized a turbulent Schmidt 
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Figure 17. Vertical Centerline Profiles at X/D =120; Case 3, Sc t = 0.2. 



Figure 18. Fuel Mixing Efficiency Distribution; Case 3, Sc t - 0.2. 



number of 0.2. 


Hydrogen contours presented in Fig. 16 illustrate fuel penetration, decay of the hydrogen 
core, and details of the upstream flow separation region. The 0.10 contour extends only about 
1.1D upstream of the jet center, and the .005 contour extends about 3. 2D upstream. The 
hydrogen core penetrates about 2.5D, and the 0.40 H 2 concentration persists only 5 jet diame- 
ters before dispersing. 

Vertical jet centerline profiles at the 120D station are compared in Fig. 17 to experimen- 
tally measured values. In addition to the high fuel concentration at the wall for the computa- 
tional results, the major differences arise from the variation in experimentally measured static 
pressure between 2 and 5 diameters from the wall. Because of the relative uniformity of the 
pressure for the other hydrogen cases, the pressure variations in the experimental data here are 
suspect. Notice the reasonable although slightly overpredicted maximum H 2 concentration and 
the extremely accurate agreement in the penetration of the .005 contour. 

Computational mixing efficiency, defined in the nonstandard way discussed above 
(Eqn. 1 and 2) for Case 3 is compared with experimental results in Fig. 18. Whereas the 
correctly defined mixing efficiency is 100%, this nonstandard mixing efficiency actually allows 
quantitative comparison of two fully mixed results. Computational results using Sc, = 0.2 
underpredict the experimental mixing efficiency by 4% and over predict the adjusted mixing 
efficiency by 8%, consistent with predictions for Case 1 which had the thin boundary layer. 

Case 4 

Only one computational solution was obtained for the 30° multiple hydrogen jets. The 
solution utilized a turbulent Schmidt number of 0.2, and results from that case are illustrated 
in Figures 19-21. 

Hydrogen contours presented in Fig. 19 illustrate fuel penetration, decay of the hydrogen 
core, and details of the upstream flow separation region. The 0.10 contour extends about 2.0D 
upstream of the jet center, and the .005 contour extends about 2.5D upstream. The core 
penetrates about 2.0D, and the 0.40 H 2 concentration persists only 5.4 jet diameters before 
dispersing. 

Vertical jet centerline profiles at the 120D station are compared in Fig. 20 to experimen- 
tally measured values. As with the cases from Ref. 4, the computed pressure underpredicts 
experimental values at the 120D station. Small differences in predicted and measured tem- 
perature and velocity are believed associated with boundary layer energization by the down- 
stream angled fuel. Note that at Z/D = 2 the velocity is higher than the freestream value. 
Peak hydrogen mass fraction occurs at the wall for the computed flow, compared to about 3D 
above the wall for the data. The computed a max is about 33% higher than experimental meas- 
urements. Based on results for Cases 1-3, this difference suggests that the predicted mixing 
efficiency for Sc, = 0.2 is about 6% lower than the experimental values. 

Longitudinal jet centerline wall pressure distribution in the nearfield is illustrated in Fig. 
21. These results show the extent of the flow disturbance upstream of the injector and a very 
small low pressure region downstream of the injector. The pressure rise extends about ID 
farther upstream than the .005 H 2 mass fraction contour. Unlike the Case 2 results, the down- 
stream pressure increases following a low pressure region. This rise in pressure is believed to 
be associated with the multiple jet interacting bow shocks. Data at 3D on the jet centerline is 
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Figure 20. Vertical Centerline Profiles at X/D =120; Case 4, Sc t = 0.2. 



Figure 21. Longitudinal Centerline Wall Pressure; Case 4, Sc t = 0.2. 
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in good agreement with prediction, however the data between jets is not. 

Case 5 

Only one computational solution was obtained for each of Cases 5 through 7. These 
solutions used Sc t = 0.5, because that value provided excellent agreement with the experimen- 
tal data. Comparisons of computational and experimental results for Case 5, 30° matched 
static pressure helium injection, are presented in Figures 22-24. 

Helium contours presented in Fig. 22 illustrate fuel penetration, decay of the core and 
details of the upstream flow separation. The .005 contour extends upstream only about 0.5D 
from the edge of the highly elongated angled injector. The He core penetrates about ID off 
the wall, and the 0.40 contour extends downstream about 8D. Basically, the helium does not 
penetrate out of the boundary layer because of the low jet-to-freestream dynamic pressure 
ratio. Case 5 had a q of 0.6 versus 1.0 for the the 30° injection of hydrogen in Case 4. 

Vertical jet centerline mole fraction profiles at 20, 40 and 90 jet diameters are compared 
with experimentally measured values in Fig. 23. The nearfield core penetrates about 1.5D off 
the plate, but the outer edge of the mixing region remains within 4D, about half as high as 
data. In the farfield, a weak peak remains about 2D off the plate, unlike the H 2 cases. Peak 
concentration are in good agreement with data throughout. 

Figure 24 presents longitudinal distributions of maximum helium mass fraction and pene- 
tration. Excellent agreement in the peak He concentration is illustrated in Fig. 24a, which was 
shown by the hydrogen cases to be a good indication that mixing efficiency is accurately 
predicted. Predicted helium penetration is quite low compared to data. Fig. 24b illustrates 
computed and experimental measurements of He penetration to the outer edge of the jet {x He 
= 0.) and the to the location of peak concentration. Measured penetration of the edge is twice 
the computed value for the nearfield and about 25% higher than computed for the 90D station. 
However, the peak penetration is accurately modeled in the nearfield and only about 30% 
underpredicted at the 90D station. 

Case 6 

Comparisons of computational and experimental results for Case 6, 30° underexpanded 
helium injection, are presented in Figures 25-27. This solution also utilized a Sc t of 0.5. 

Helium contours presented in Fig. 25 illustrate fuel penetration, decay of the core and 
details of the upstream flow separation. The .005 contour extends upstream about ID from 
the edge of the highly elongated angled injector. The He core penetrates about 2.5D off the 
wall, and the 40% contour extends downstream past the end of the elliptic solution (>10D). 
Comparing Fig. 25 with Fig. 22, the underexpanded jet has much higher penetration, con- 
sistent with the fivefold increase in helium pressure and mass flow rate. 

Vertical jet centerline mole fraction profiles at 20, 40 and 90 diameters are compared 
with experimentally measured values in Fig. 26. This high pressure jet penetrates about 4D 
off the plate in the nearfield, and the computational profile maintains a peak concentration 6D 
off the plate at the 90D station. Shape and magnitude of the x He contours are in excellent 
agreement with experimental measurements at all three stations. The only noticible 
discrepancy is near the wall, where the experimental profile drops off more rapidly than the 
computational results. 
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Figure 22. Longitudinal Helium Mass Fraction Contours; Case 5,9 = 30° , q = 0.6. 



a) X/D = 20 b) X/D = 40 c ) X/D = 90 

Figure 23. Vertical Centerline Helium Mole Fraction Profiles; Case 5, 9 = 30°, q = 0.6. 
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Figure 24. Longitudinal Distributions; Case 5, 9 = 30°, q = 0.6. 
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Figure 25. Longitudinal Helium Mass Fraction Contours; Case 6, 0 = 30°, q = 3.1. 
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Figure 26. Vertical Centerline Helium Mole Fraction Profiles; Case 6, 0 = 30°, q =3.1. 
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Figure 27. Longitudinal Distributions; Case 6, 0 = SO 0 , q = 3.1. 
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Figure 27 presents longitudinal distribution of maximum helium mass fraction a max and 
penetration P/D. The excellent agreement in Fig. 27a indicates that mixing efficiency is accu- 
rately predicted. Computed helium penetration, both core and outer edge of mixing region, are 
compared to experimental results in Fig. 27b and show excellent agreement. The largest 
discrepancy is the 15% underpredicted penetration of the edge at X/D = 20. 

Case 7 

Comparisons of computational ( Sc t = 0.5) and experimental results for Case 7, 15° 
matched static pressure injection, are presented in Figures 28-30. 

Helium contours presented in Fig. 28 illustrate fuel penetration, decay of the core and 
details of the upstream flow separation. The .005 contour extends upstream about 0.5D from 
the edge of the highly elongated angled injector. The He core penetrates about 0.5D off the 
wall, and the 40% contour extends downstream about 10D. 

The low pressure, low angle injector produces a mixing region characterized by the 
highest He concentration and core on the wall by the first survey station, X/D =20, as illus- 
trated by the vertical profiles in Fig. 29. The computational results indicate that a core exists 
off the wall in the nearfield, but returns to the wall by the 20D station. Shape and magnitude 
of the predicted mole fraction contours are in good agreement with experimental contours. 
However, the underprediction of the peak x He by about 40% at the 90D station suggests that a 
slightly lower value of Sc t may be appropriate for this case. However, this is inconclusive, 
since the experimental data for this case showed an unexplained asymmetry. 

Figure 30 presents longitudinal distribution of maximum helium mass fraction a max and 
penetration P/D. Fig. 30a illustrates good agreement between calculated and measured peak 
helium mass fraction, suggesting that the computational results accurately predict mixing 
efficiency. Predicted helium penetration to the outer edge of the mixing region is compared to 
the experimental data in Fig. 30b, illustrating slight underprediction at 20 and 40 X/D, but 
excellent agreement at the 90D station. 

DISCUSSION OF RESULTS 

Qualitatively, the computations are good in that computational trends are consistent with 
experimental observations. The results show that increasing 8/D increases penetration, both of 
the outer edge of the mixing region and of the point of maximum concentration, which was a 
conclusion of Ref. 4. Pressure plots (not included) reveal bow shocks which are naturally 
stronger for the normal injection and underexpanded injection. Velocity vector plots (also not 
included) reveal the horseshoe vortex system characteristic of injection flow fields. These and 
other observations and details of the flow solutions are beyond the scope of this presentation, 
but will be addressed in the future. 

Accuracy 

Table 4 summarizes observed differences between computed and experimental measure- 
ments for jet penetration (P/D and Z m [D), spreading (W/D and W B fD), maximum H 2 concen- 
tration (a max ), and mixing efficiency, v \ m . 

Where available, solutions using Sc t = 0.5 as recommended above are included in the 
table. Because of the recognized errors in the hydrogen data, the Sc t = 0.2 solutions are used 


662 




Figure 28. Longitudinal Helium Mass Fraction Contours; Case 7, 0 = 15°, q = 0.27. 
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Figure 29. Vertical Centerline Helium Mole Fraction Profiles; Case 7, 0 = 15°, q = 0.27. 



a) Maximum Helium Mass Fraction b ) Helium Penetration 

Figure 30. Longitudinal Distributions; Case 7, 0 = 15°, q = 0.27. 
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when comparing peak hydrogen mass fraction and mixing efficiency. Mean deviation is sum- 
marized in the last column. 

The codes consistently underpredicted penetration of the maximum concentration (Z m /D) 
in the farfield, generally leaving the peak along the wall. Penetration to the outer edge of the 
mixing layer is generally predicted within ±1.6 diameters, or to within about ±7%. 

Lateral spreading is slightly overpredicted in the freestream and generally significantly 
overpredicted in the boundary layer. Comparisons available are insufficient to provide reliable 
estimates, but the range is about ID for W and -2D for W B . 

Peak injectant concentration in the farfield is over or underpredicted by about ±17%, 
while the unadjusted fuel mixing efficiency has an average deviation of 5% (universally under- 
predicted). Confidence in this mixing efficiency accuracy will be significantly enhanced by 
additional data and integration of the available helium data. 

Table 4. Accuracy Summary 


Case no. 



1 

2 

3 

4 

5 

6 

7 

M.D. 

X/D 

120 

120 

120 

120 

90 

90 

90 


Sc, 

0.2 

0.5 

0.2 

0.2 

0.5 

0.5 

0.5 


P/D 

2.0 

4.0 

0.0 

3.0 

-1.8 

-0.5 

0 

1.6 

ZJ D 

-2.5 

-3.0 

-2.5 

-2.5 

-1.5 

0.2 

0 

1.7 

W/D 

1.25 

.70 

— 

— 

— 

— 

-- 

.98 

W B / D 

4.0 

0.2 

__ 

— 

— 

-- 

— 

2.1 

Cf 

^max 

-14% 

0%* 

20% 

30% 

10% 

12% 

-30% 

17% 


-.09 

-.05* 

-.03 

— 

— 

— 

— 

.05 


*Sc, =0.2 


Recommendations 

There are obvious differences between the hydrogen and helium cases in the application 
of the code to achieve acceptable code vs. experiment agreement. For the H 2 cases, better 
agreement was obtained using Sc, = 0.2, whereas for the He cases, better agreement was 
obtained with Sc, = 0.5. Because of the superior experimental measurement capability of the 
recent helium tests reported in Ref. 30 and documented uncertainties of the older hydrogen 
data[3,4 and 7], it is prudent to recommend solution techniques used for Cases 5-7. There- 
fore, when using a fine grid and the Baldwin-Lomax turbulence model, results indicate that 
mixing of transverse injectors with injection angles, between 15 and 90 degrees should be 
modeled using Sc,~ 0.5. However, when performing a complete combustor analysis requiring 
a coarse grid with multiple fuel injectors, a more realistic prediction of fuel mixing may be 
achieved using Sc, - 1.0[28]. It would also be advisable to compare grid density used for 
such solutions to the injection cases (preferably He) presented in this paper. 

Fuel mixing efficiency is relatively insensitive to maximum fuel concentration (Cases 1 
and 2). If the decay of peak fuel concentration is accurately predicted, then the mixing 
efficiency will be accurate, but not vice versa. 
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Turbulent diffusivity (Schmidt number) has been shown to vary across a boundary 
layer[37] and that trend appears in Table 4, Lateral spreading within and spreading out of a 
boundary layer are not correcdy predicted when using a single value of Sc t over the entire 
flow. This observation suggests that improved turbulence modeling to account for the nonuni- 
form turbulent diffusivity is in order. 


CONCLUSIONS 

The SPARK family of Navier Stokes codes were validated for fuel mixing in a super- 
sonic flow, using experimental data for discrete circular injectors angled at 15° (downstream) 
to 90° (normal) for air flow boundary layer thicknesses ranging from 1.25 to 6.25 injector 
diameters. Details of the physical model, solution methodology, boundary conditions and tur- 
bulence model are presented in sufficient detail such that recreation of the solutions is possi- 
ble. Results presented demonstrate that "grid converged" solutions for one injector were not 
attained but that mixing efficiency sensitivity is acceptable. Mixing efficiency is also rela- 
tively insensitive to peak fuel concentrations. Recommendations are presented for modeling 
the fuel injector problem. These recommendations are different depending on the grid density 
around each injector. Coarse grids require damped turbulence (achieved with a turbulent 
Schmidt number of unity). This type of grid is representative of large scale combustor solu- 
tions with many fuel injectors. Fine grids require a reduced turbulent Schmidt number on the 
order of 0.5, which are in agreement with experimentally observed values for similar flows. 
The turbulent Schmidt number is not constant across the entire flow domain; a distribution 
may enhance fuel injection modeling. 
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